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ABSTRACT 

We analyze the statistical properties of the turbulent velocity field in the deflagration model for Type la 
supernovae. In particular, we consider the question of whether turbulence is isotropic and consistent with 
the Kolmogorov theory at small length scales. Using numerical data from a high-resolution simulation of a 
thermonuclear supernova explosion (Ropke et al., 2007), spectra of the turbulence energy and velocity struc- 
ture functions are computed. We show that the turbulent velocity field is isotropic at small length scales and 
follows a scaling law that is consistent with the Kolmogorov theory until most of the nuclear fuel is burned. 
At length scales greater than a certain characteristic scale that agrees with the prediction of Niemeyer and 
Woosley (1997), turbulence becomes anisotropic. Here, the radial velocity fluctuations follow the scaling law 
of the Rayleigh-Taylor instability, whereas the angular component still obeys Kolmogorov scaling. In the late 
phase of the explosion, this characteristic scale drops below the numerical resolution of the simulation. The 
analysis confirms that a subgrid-scale model for the unresolved turbulence energy is required for the consistent 
calculation of the flame speed in deflagration models of Type la supernovae, and that the assumption of isotropy 
on these scales is appropriate. 

Subject headings: hydrodynamics — instabilities — methods: statistical — turbulence — Supernovae: general 



1. INTRODUCTION 

The mechanism of thermonuclear explosions of white 
dwarf (WD) stars, giving rise to Type la supernovae (SNe la) 
is still not understood in full detail. For three reasons, a key to 
this problem is the understanding of turbulent thermonuclear 
combustion in the deflagration phase — the phase of subsonic 
flame propagation which commences the explosion process. 
First, it is required for the correct modeling of the flame prop- 
agation in the deflagration model of SNe la (see |Hilleb randt 
& Niemeyer 2000 for a review of SN la explosion scenar- 
ios). Second, in alternative models, burning starts out in the 
deflagration mode, and this is an essential ingredient to the 
overall explosion process. Third, turbulence in the deflagra- 
tion phase sets the conditions for a possible deflagration-to- 
detona tion transition (DDT) in the delayed detonation sce- 
nario (Ropke 2007 Woosley 2007). The necessary insight 
into the details of the turbulent combustion process, however, 
is hampered by the fact that full-star simulations of thermonu- 
clear supernova explosions cannot resolve the structure of the 
deflagration flame. At the large scales accessible to simula- 
tions, the flame propagation is determined by flame instabil- 
ities and turbulence. These effects significantly boost the ef- 
fective burning speed and, in the pure deflagration model of 
SNe la, lead to the flame acceleration required to explode the 
WD ( |Reinecke et al.|2002| |Gamezo et al.|2003| >. In order to 
describe the flame propagation in such simulations, the in- 
teraction of the flame with turbulent velocity fluctuations has 
to be modeled. These modeling approaches yield an effective 
flame propagation speed on the numerically resolved scales — 
the so-called turbulent burning speed. 

The problem of calculating this turbulent flame propagation 
speed in three-dimensional simulations of thermonuclear su- 



pernova explosions 1 has been the subject of a lively debate. 
One school of thought holds the view that the effective prop- 
agation speed in the flamelet regime would naturally be given 

by the velocity scale vrt(^) oc (gegl) 1 ^ 2 associated with the 
Rayleigh-Taylor (RT) instability induced by buo yancy in the 
gravitational field g e s for any length scale £ (Sharp 1984]!. As 
a subgrid scale model (I = A, where A is the numerical cutoff 
length), this scaling relation is easily implemented and ap- 
pears to be motivated by the basic physics of thermonuclear 
combustion in Type la supernovae (Gamezo et al. 20031. In 
opposition to this vie w, |Niemeyer & Hillebrandt| ( |1995f and 
|Niemeyer & Kerstein| ( |1997[ ) argued that inevitably turbulent 
velocity fluctuations v (£) are dominated by the turbulent cas- 
cade at length scales I small compared to the scale of en- 
ergy injection by the RT instability an d, hence, should follow 
the Kolmogorov scaling v'(£) oc . Niemeyer & Woosley 



( 1997 1 estimated the transition length £ K / RX between the RT- 
dominated length scales (the "large scales" I > ^k/rt) ar, d the 
regime of the turbulent cascade (the "small scales" £ < ^k/rt) 
to be of the order 10km. Since the cutoff length A in con- 
temporary numerical simulation is comparable to or less than 
£k/rt> it follows that a subgrid scale (SGS) model for the con- 
sistent calculation of the turbulent flame speed has to be based 
on the turbulence energy associated with the length scale A, 
which is determined by the the dynamics of the turbulent cas- 
cade. Such an SGS model was proposed by Niemeyer & 
Hillebrandt| ( |1995 1 and further developed by Schmidt et al. 
( |2006[ ). iTTs not clear a priori that the turbulence to be cap- 
tured by the SGS ansatz is of Kolmogorov-type and therefore 
the approach of Schmidt et al. (2006) is not based on this as- 
sumption. However, Kolmogorov scaling is an obvious possi- 



For a review, see Ropke & Schmidt 2008 . 
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bility that has to be considered. 

The best way to gain insight into the properties of turbu- 
lence in the deflagration stage is to analyze it directly in high- 
resolution simulations of the deflagration model. Here the ef- 
fects of gravity and spherical expansion of the background 
are naturally accounted for. In this paper, we present such 
an a posteriori analysis of turbulence based on data from a 
recent numerical simulation (Ropke et al. 2007) which was 
carried out on a very large grid (1024 J cells). The character- 
istics of this deflagration model as well as derived synthetic 
observables are in reasonable agreement with the observa- 
tions of dim mer (but still normal) observed SNe la (Ropke 
et al.|2007] i. Therefore, our analysis is based on data from a 
model that is expected to give a realistic picture of turbulence 
in SNe la. The failure of the pure deflagration scenario to re- 
produce the brighter end of the SN la sample does not limit 
the significance of our results as in all alternative scenarios 
currently under discussion a similar deflagration phase initi- 
ates the explosion and sets the stage for the later evolution. 

In the simulation analyzed here, the co-moving grid tech- 
nique introduced by Ropke ( 2005 ) allowed for a very small 
initial cutoff length A ~ 1km in the inner regions of the ex- 
ploding white dwarf. Although A was gradually increased 
in the course of the simulation, it was possible to investi- 
gate the behaviour of turbulent velocity fluctuations at length 
scales I ~ 10km by means of computing kinetic energy spec- 
trum functions of the velocity field in subdomains selected by 
an appropriate window function. The results indicated Kol- 
mogorov scaling (Ropke et al. 2007). In this article, we re- 
fine this analysis by decomposing the velocity field into ra- 
dial and angular components. In addition, we subtract the 
spherically averaged radial velocity in order to separate the 
velocity fluctuations from the mean radial expansion of the 
white dwarf. For both components of the fluctuating veloc- 
ity field, a Fourier analysis is carried out to compute kinetic 
energy spectra and to investigate possible anisotropies. Since 
Fourier transforms are distorted by the non-periodic boundary 
of the computational domain, we calculate structure functions 
of the fluctuating velocity field to obtain reliable estimates of 
the scaling properties for the whole dynamical range of the 
simulation. 

The methodology of our analysis is explained in detail in 
the following section. As will be shown in Section|3] RT scal- 
ing is found for the radial fluctuating velocity component at 
length scales greater than £k/rt, whereas Kolmogorov scal- 
ing applies for smaller length scales. In close agreement with 
the prediction by Niemey er & Woosley| ( |1997] l, the numeri- 
cally determined value of ^k/rt is about 15 km after the on- 
set of the explosion. The angular velocity component, on the 
other hand, closely follows the Kolmogorov scaling law at 
all length scales. Furthermore, the velocity fluctuations are 
nearly isotropic at length scales smaller than £k/rt- In the l ast 
section, we discuss possible caveats of our analysis and com- 
ment on the implications for numerical simulations of Type la 
supernovae. 



2. ANALYSIS OF THE TURBULENT VELOCITY FIELD 

For the statistical analysis of turbulence, we have to take 
into account that the velocity field is a superposition of turbu- 
lent velocity fluctuations and the bulk expansion of the white 
dwarf. To estimate the bulk expansion, we average the radial 
component of the velocity field over spherical shells of dis- 



crete radii r, : 



1 Nl 

v(n) = — ■ 2jv(rf/) • e r (r ,;), 



(1) 



where rf = iA 2 (t) is an integer multiple of the squared size 
A(f) of the grid cells at time t. The sum is over all cells in the 
cubic grid that are located at the distance r, from the center, 
and Nj is the corresponding number of cells. The unit vector 
in radial direction at the position r, j of the j-th cell in the i- 
th shell is denoted by e,-(r,7). Using this estimate, we subtract 
the spherically averaged component v(r,) from the original ve- 
locity field. In the following, it is understood that the symbol 
v refers to the fluctuating part of the velocity field. 

Since the RT-Instability evolves in the direction of gravity, 
we perform all computations with velocity components paral- 
lel and perpendicular to the gravitational field. These compo- 
nents correspond to the radial and angular directions, because 
the gravitational field is assumed to be spherically symmetric 
in the simulation. Thus, we define vii := v r e r , where e r = r/r 
and v r = v ■ e f , and vj_ := v — vy . The corresponding energy 
spectrum functions are obtained by integrating the kinetic en- 
ergy per unit mass over spheres of radius k in Fourier space: 



£||(fc)=i/d^ 2 | V ||(k)| 2 , 

E ± (k)=tfdn k k 2 \\ ± (k)\ 2 , 



(2a) 
(2b) 

where Vii(k) and vj_(k) are the Fourier transforms of the lon- 
gitudinal and transversal velocity components, respectively. 
For developed turbulence, it is expected that the the energy 
spectrum functions follow power laws, E(k) oc k~@, in the in- 
ertial subrange of wave numbers. A disadvantage of Fourier 
transforms is that the contributions from small wave numbers 
are distorted by the non-periodic boundaries of the computa- 
tional domain. For this reason, we apply Gaussian window 
functions to the data sets as described in Ropke et al. (2007). 
Since the data windowing corresponds to a high-pass filter in 
Fourier space, the range of the energy spectra is constrained 
to higher wave numbers. 

In contrast to the energy spectra, structure functions are 
two-point velocity correlation functions computed in position 
space. While the Fourier transforms are performed in Carte- 
sian coordinate systems, we use spherical coordinates for the 
computation of structure functions, which is convenient to de- 
fine directions parallel and perpendicular to gravity. More- 
over, the computation can be constrained to the interiors of 
spheres containing a certain amount of burned matter. We de- 
fine the radial velocity increment by the difference of vii , i. e., 
the velocity component in the direction of gravity, at two dif- 
ferent positions rj and r 2 , 



<5v ra 



rati ' 



i(r 2 )-V||(ri) 



(3a) 



The angular velocity increment is defined by the difference 
of the velocities projected in the directions perpendicular to 
gravity, i e., 



£v an g = v_L(r 2 )-v_L(ri) 



(3b) 



Note that these increments do not correspond to longitu- 
dinal and transversal velocity increments, because the veloc- 
ity components are generally not parallel or perpendicular to 
the spatial separation r 2 -ri. However, we think that the 
above definitions of velocity increments are better suited to 
the physics of RT-driven turbulence in thermonuclear super- 
novae. Our proposition is corroborated by the the scaling 
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properties of the structure functions that will be presented in 
the following Section. 

The radial and angular structure functions of order p are 
defined by the averages of the radial and angular velocity in- 
crements to the power p, respectively: 



VadW=(I^V rad H 



(4a) 
(4b) 



where the length scale £ := |ra — ri|. There is a large range 
of length scales which encompasses the turbulent interior of 
the exploding WD. For fully developed turbulence the struc- 
ture functions are given by power laws Sp jIa d(i) oc ^'"" d and 
S p ^ ng (£) oc v* p - ms , where Cp,rad and Cp.ang are characteristic scal- 
ing exponents. If we further assume isotropy, C^.rad — C 



Kol- 



and C, p = p/3 according to the theoretical analysis by 
mogorov| (1941) . In particular, it follows that the turbulent 
velocity fluctuation v'(£) oc £ 1 ^ . For the turbulent flow driven 
by the RT insta bility, on the other hand, v'{£) oc t 1 / 1 ( |Davies 
& Taylor 1950) corresponding to Q p = p/2. 

For the numerical computation of the structure functions, 
one has to take a sufficient large number of sample points 
that are distributed with uniform probability within a spher- 
ical region of prescribed radius in order to achieve converged 
statistics. This was achieved by a Monte-Carlo-type algo- 
rithm, where the total number of sample points was varied 
and, thereby, convergence was established. 

In order to analyze the isotropy of the velocity field close 
to the flame, we performed calculations in small boxes inter- 
sected by the flame. The algorithm is based on the analysis 
performed by Zing ale et al.| p005). To simplify the calcula- 
tion, the box is placed along a coordinate axis, in our case the 
z-axis, which defines the local direction of gravity. For each 
cell within the box, the velocity difference S\ between the 
local velocity and the velocity at the center of the box is cal- 
culated. Then projected contours of the velocity differences 
in Fourier space can be constructed, by integrating (5v(k) over 



circles of radius k p = ^Jk 2 + k 2 . in planes perpendicular to the 

z-component of k. This procedure was applied for several po- 
sitions of the box center corresponding to different fractions 
of burned matter in the box. 

3. RESULTS 

3.1. Energy spectra 

We plot the energy spectra as functions of the normalized 
wave number k n = 512A ( )(f)A:/7r for 16 < k n < 512, where 
Ao(f) is the size of the cells in the uniform part of the grid at 
time t. We have A (f) = 2.93 and 14.69km at times t = 0.5 and 
1.0 seconds, respectively. Note that wave numbers k n < 32 
are obscured by data windowing (see Section [2]). The com- 
puted energy spectrum functions at 0.5 seconds are shown as 
black curves in Fig. [T] (a). The thin gray lines corresponds 
to power-law fits and the dashed line indicates the expected 
spectrum according to the Kolmogorov theory with an expo- 
nent (3 = 5/3. For the longitudinal spectrum function, we find 
j3 ~ 1 .64 for high wave numbers, which is in good agreement 
with the Kolmogorov theory. The longitudinal spectrum pos- 
sibly stiffens toward lower wave numbers, but the accuracy 
of the computed spectra does not allow for a conclusive re- 
sult. We exclude the lower part of wavenumbers from the 
fit, because the slope becomes steeper than a Kolmogorov 
spectrum for k a < 64. The transversal spectrum is slightly 



steeper (f3 ~ 1-71), but close to a Kolmogorov spectrum for 
the whole range of wave numbers. These results were corrob- 
orated by the computation of second-order structure functions 
with much higher accuracy (see section |3.2[ ). We also note 
that the Kolmogorov scali ng for higher w ave numbers agrees 
with the findings of Ropke et al.l (|2007[l, where turbulence 



energy spectra were computed without splitting the velocity 
field. This is in accordance with the expectation that the ra- 
dial expansion will mostly affect low wave numbers modes. 
At t = 1 .0 seconds, on the other hand, the longitudinal spec- 
trum has an exponent (3 « 1 .97 over the entire range of wave 
numbers. In contrast, the transverse spectrum is much shal- 
lower. The exponent (3 = 2 corresponds to RT scaling, because 
E(k) oc k~ 2 implies 8v(£) oc £}l 2 . Consequently, it appears that 
the velocity component in the direction of gravity is domi- 
nated by the RT instability even at the smallest numerically 
resolved scales in the late phase of the explosion, while Kol- 
mogorov scaling is found for the velocity component perpen- 
dicular to gravity. 

3.2. Velocity structure functions 

For the computation of the structure functions following 
equation pa] ) and ( HE} , we chose a spherical region in which 
90% of the material was burned. This choice of the re- 
gion resulted from the requirement of encompassing the bulk 
of turbulence at a given time, while excluding the outer, 
non-turbulent regions of the white dwarf. Fig. [2] shows 
double-logarithmic plots of the radial (solid curve) and an- 
gular (dashed curve) structure functions of second order at 
different times. In all cases up to t = 0.7 seconds, a scal- 
ing exponent £2 close to 2/3 is found for the range of length 
scales I < 10km. Moreover, S p ^(£) ~ S p ^ ng (£), which indi- 
cates isotropy at small length scales. For larger length scales, 
on the other hand, the radial structure functions S p . mt i(£) obey 
a scaling law with an exponent Ci.md « 1, whereas S Ptimg (£) 
still follows Kolmogorov scaling. From C2,md ~ 1, it follows 
that v'(£) oc £ l l 2 . As outlined in section^ this corresponds to 
RT scaling. For this reason, the change of slope of the radial 
structure function indicates a transition from the inertial-range 
turbulence cascade to the regime of RT instabilities at a length 
scale ^k/rt approaching w 14 km in the course of the explo 
sion. Thus, our analysis confirms the estimate by Niemeyer 
|& Woosley| ( |1997[ l. As a result of the overall expansion of the 
co-moving grid, the transition length drops below the numer- 
ical resolution of the simulation after ~ 0.7 to 1.0 seconds. 
The time evolution of ^k/rt is further illustrated in Fig. [3j 
where the second order structure functions are plotted for an 
data sets, for which £k/rt is numerically resolved. While in 
Fig. [2]the range of length scales is adjusted to the size of the 
co-moving grid, a fixed range of length scales is used in Fig. [3] 
Additionally, the typical mass density in the vicinity of the 
flame front is specified for each instant of time. In agreement 
with Fig. 1 of |Niemeyer & Woo sley ( 1997), Fig. 2 and 3 show 
that £k/rt becomes smaller with decreasing density (and ad- 
vancing time). Remarkably, it appears that ^k/rt approaches 
an asymptotic value, but we cannot investigate this behavior 
for t > 0.8 seconds. To obtain more properties of £k/rt ad- 
ditional high-resolved numerical simulations are needed, in 
which £k/rt can be tracked for a longer time. 

In Figure EJ the structure functions up to the sixth order are 
plotted for t = 0.5, 0.6 and 0.7 seconds. As one can see from 
the scaling exponents listed in Figure [3] C3 is close to unity 
for the angular structure functions and for the radial structure 
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FIG. 1 . — Longitudinal (1) and transversal (2) energy spectrum function (black curve) at t = 0.5 (a) and t = 1.0 (b) seconds in comparison with the Kolmogorov 
energy spectrum (dashed line). 
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FIG. 2. — Radial (solid curve) and angular (dashed curve) 2nd-order structure function with the determined transition length scale. 



functions in the subrange I < Ik /rt- This res ult is consistent 
with the Kolmogorov theory (see Frisch|1995|l. For £ > £k/rt> 
on the other hand, (3 ,3d rj 1.5. KemarkaDly, it appears that 
for all p < 6 the slopes of the radial structure functions are 
steeper by a factor of 1.5 at length scales greater than ^k/rt- 
This suggests that the ratio Z p := Cp/d is approximately equal 
for S p , ad (t < £ K / RT ) and S p ^ A (i > £ K / RT ). Indeed, the rel- 
ative scaling exponents of the radial structure functions in 
both subranges are quite close (see Table [TJ. There are the- 
oretical arguments (Dubrulle 1994) as well as numerical in- 
vestigations ptenzi et al.||1993| l in support of the fundamen- 



tal significance of relative scaling exponents. But it has not 
been noticed before that RT-driven velocity fluctuations ex- 
hibit statistical properties that are equivalent to the proper- 
ties of isotropic, inertial-range turbulence in terms of relative 
scalings. It is known that the calculation of higher-order ex- 
ponents becomes increasingly uncertain due to sampling er- 
rors ([Frisch 1995). For this reason, we do not consider the 
relatively high discrepancies between the results for p > 5 to 
be significant. Comparing the relative scaling exponents of 
the radial and the angular structure functions in the subrange 
t > £k/ki, on the other hand, we find very good agreement. 
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TABLE 1 

Relative scaling exponents Z p at three different instants of time 







fM 


p = l 


p = 2 


p = 3 


p = A 


p = 5 


p = 6 






0.5 


0.344± 0.016 


0.676 ±0.027 


1.0 ±0.037 


1.309 ±0.050 


1.597 ±0.058 


1.860 ±0.068 


Z Pl rad for < 


< <?K/RT 


0.6 


0.341 ±0.016 


0.676 ±0.030 


1.0 ±0.046 


1.304 ±0.057 


1.593 ±0.071 


1.850 ±0.082 




0.7 


0.347 ±0.019 


0.676 ±0.029 


1.0 ±0.038 


1.302 ±0.052 


1.584 ±0.060 


1.842 ±0.074 






0.5 


0.339 ±0.009 


0.674 ±0.015 


1.0 ±0.021 


1.315 ±0.027 


1.615 ±0.033 


1.900 ±0.038 


Z,,,md for t 


> ^K/RT 
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0.338 ±0.011 
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1.0 ±0.027 


1.316±0.033 
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1.896 ±0.050 
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0.336 ±0.008 


0.670 ±0.016 


1.0 ±0.024 


1.323 ±0.034 


1.636 ±0.042 


1.920 ±0.050 
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FIG. 3 . — 2nd order radial structure functions with the transition length (big 
dots) and corresponding density on the flame front. 

We will concentrate on these scaling exponents in the follow- 
ing. 

3.3. Local isotropy of the velocity fluctuations 

The contour lines of the Fourier-transformed velocity field 
at 0.5 seconds inside a box as described in Section |2] are 
plotted in Fig [5] The box is positioned such that 50% of 
the enclosed matter is burned. The thick dashed line indi- 
cate the wave number 27r/£ K / RT corresponding to the transi- 
tion length obtained from the calculation of the radial struc- 
ture functions. For smaller wave numbers (i. e., I > £ RT ), the 
contours of 5v(k) are clearly anisotropic. One can see that 
a given velocity difference in the z-direction spans a smaller 
range of wave numbers than in the directions perpendicular 
to the z-axis. This corresponds to a steeper slope of the ve- 
locity fluctuations in the radial direction, which is approxi- 
mately given by the z-direction. Toward higher wave num- 
bers (£ < ^ RT ), the anisotropy of the contours decreases, but 
we do not find perfect isotropy. This might indicate residual 
small-scale anisotropy in the vicinity of the flame front, which 
can be caused by the intermittency of turbulence. However, 
it could also be a spurious effect due to the misalignment be- 
tween the z-axis and the radial direction at off-center positions 
within the box. 

4. CONCLUSION 

We investigated the scaling properties of turbulence in a 
high-resolution simulation of a Type la supernova based on 
the pure deflagration model (Rtipke et al .|2007 1. Both energy 



spectrum functions and structure functions of second order 
were computed. The results of this study are as follows: 

1. The velocity fluctuations in the radial direction, i. e., 
the direction of gravity, follow Kolmogorov scaling at 
length scales smaller than a certain transition length 
£k/rt ~ 10km. Only at length scales greater than ^k/rt> 
the radial velocity fluctuations are dominated by the 
scaling law of the Raylei gh-Taylor instability. T his be- 

1997]). 



haviour was predicted by Niemeyer & Woosley 



2. The velocity fluctuations in angular directions, i. e., per- 
pendicular to gravity, obey Kolmogorov scaling over 
the entire range of numerically resolved length scales. 

3. For £ < 

^k/rtj the magnitudes of the radial and angu- 
lar velocity fluctuations are nearly equal. For this rea- 
son, small-scale turbulence appears to be statistically 
isotropic. 

4. Fourier analysis of the velocity fluctuations in a small 
region near the flame surface allows for slight residual 
anisotropics at the smallest resolved scales. 

As regards the interpretation of our results, a possible cause 
for concern is that Kolmogorov scaling in the radial direction 
is only found for a relatively narrow range of length scales 
greater than the numerical cutoff length. It is known that 
these scales are affected by numerical dissipation and, par- 
ticularly, by the bottleneck effect (Schmidt et al. 2006). Thus, 
the flattening of the turbulence energy spectrum might be ar- 
tificial. However, in this case, no significant flattening should 
be observed for the corresponding structure functions which 
are much less affected by the bottleneck effect. The scal- 
ing laws implied by the radial energy spectra and structure 
functions are fully consistent at small length scales and the 
transversal spectra show no flattening at all. In consequence, 
we are confident that the scaling laws are genuine. A possi- 
ble explanation for the absence of the bottleneck effect is that 
turbulence does not reach a statistically stationary state in a 
supernova explosion. 

Yet another issue is that Kolmogorov scaling might be en- 
forced by the subgrid scale model that was used in the sim- 
ulation. However, other than the RT-based models used by 
|Gamezo et al.| ([2003 ), this SGS model does neither presume 
any given scaling of turbulence nor does it imprint such a scal- 
ing on the numerically resolved flow. A potential problem for 
the SGS model is the possible lack of isotropy near the flame 
surface. However, there is certainly no pronounced anisotropy 
at the smallest resolved scales if Ao(f) < £k/rt> and statistical 
isotropy is found for the bulk of turbulent regions. Only in the 
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FIG. 4. — Radial and angular structure function up to the sixth order with the corresponding scaling exponents. For the radial structure functions, the transition 
length from Kolmogorov to Rayleigh-Taylor scaling ^k/rt ' s indicated. 
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FIG. 5. — Contour plots of the Fourier-transformed velocity differences 
inside a small box intersected by the flame front at 0.5 seconds. 



late phase of the explosion, when the transition length «?k/rt 
becomes smaller than the numerical resolution and the re- 
solved small-scale turbulence definitely becomes anisotropic, 
the notion of SGS turbulence energy cannot be strictly jus- 
tified. One should note, however, that this point more or 
less coincides with the time when a deflagration-to-detonation 



transition is expected to occur ( |Gamezo et al. 2005 ; R opke &| 
Niemeyer 2007 ). Apart from that, the production of turbu- 
lence energy by unresolved buoyancy effects is heuristically 
included in the SGS model. In conclusion, the SGS turbu- 
lence energy model by Schmidt et al. ( 2006[) applies to the 



major part of the explosive burning in the deflagration phase 
of a Type la supernova explosion, but there is no regime for 
which a pure RT-scaling model holds. 

After settling the issue of turbulence scaling in the defla- 
gration phase of a Type la supernova in the present article, 
we mention that the occurrence of deflagration-to-detonation 
transitions can be constrained on the basis of the deflagration 
model. For the DDT mechanism to operate, strong turbulence 
is necessary in late phases of the burning (e.g. |Woosley|2007] ). 
Ropke (2007) found that this may indeed be realized in defla- 



gration models of SNe la with low (but not vanishing) prob- 
ability. In order to better quantify this intermittency effect, 
higher-order structure functions have to be computed. Fitting 
intermittency models to the numerically determined scaling 
exponents (as proposed by |Pan et al. 2008 ), the probability 
of strong turbulent velocity fluctuations at any instant of time 
can be estimated. This analysis will be presented in a future 
publication. 
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